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Abstract 

We study the characteristics of the outflow in disks of Be stars, based on the viscous decretion disk 
scenario. In this scenario, the matter ejected from the equatorial surface of the star drifts outward because 
of the effect of viscosity, and forms the disk. For simplicity, we adopt the a-prescription for the viscous 
stress, and assume the disk to be isothermal. Solving the resulting wind equations, we find that a transonic 
solution exists for any value of a. The sonic point is located at r > 100R for plausible values of parameters, 
where R is the stellar radius. The sonic radius is smaller for higher temperature and/or larger radiative 
force. We also find that the topology of the sonic point is nodal for a > 0.95, while it is of saddle type for 
a<0.9. We expect that the sonic point in the former case is unstable, whereas that in the latter case is 
stable. The outflow is highly subsonic in the inner part of the disk. Roughly, the outflow velocity increases 
linearly with r and the surface density decreases as r~ 2 . Interestingly, the disk is near Keplerian in the 
inner subsonic region, while it is angular momentum conserving in the outer subsonic region and in the 
supersonic region. Our results, together with the observed range of the base density for Be star disks, 
suggest that the mass loss rate in the equatorial region is at most comparable with that in the polar region. 
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1. Introduction 

Be stars are non-supergiant early-type stars with 
Balmer emission lines, whose spectral types range from 
late O- to early A-type. Extensive studies have revealed 
that a Be star has a two-component extended atmo- 
sphere, a polar region and a cool (<~ 10 4 K) equatorial 
disk. The polar region consists of a low-density, fast out- 
flow emitting UV radiation. The terminal velocity of 
the polar wind is of the order of 10 3 kms _1 . The wind 
i structure is well explained by so-called line-driven wind 
model, in which the radiative acceleration results from 
the scattering of the stellar radiation in an ensemble of 
spectral lines (Castor et al. 1975; Abott 1982). The mass 
loss rate of this region inferred from the UV lines is about 
10~ n - lO^M^yr- 1 (Snow 1981). 

On the other hand, the equatorial disk consists of a 
high-density plasma from which the optical emission lines 
and the IR excess arise. The radial velocity of the disk 
is smaller than a few kms -1 , at least within ~ 10 stellar 
radii (Hanuschik 1994, 2000; Waters, Marlborough 1994). 
The small radial velocity and the small pressure gradient 
force, as in viscous accretion disks, suggest that the equa- 
torial disk is pressure-supported and geometrically thin 
in the vertical direction and is rotationally-supported in 
the radial direction. In fact, the presence of near Keple- 
rian disks around Be stars is supported observationally as 
well as theoretically. For example, the relation between 



the size of the Ha-emitting region of Be disks resolved 
with the optical interferometer (Quirrenbach et al. 1997) 
and the separation of the double peaks of the Ha line is in 
agreement with that expected for near Keplerian disks. 
In addition, the long-term line-profile variations are well 
explained by global disk oscillations with a one-armed 
density perturbation pattern, which can be present only 
in near Keplerian disks (Okazaki 1991; Papaloizou et al. 
1992; Hanuschik et al. 1995). The mass loss rate through 
the disk is not known. The density in the disk derived 
from the IR excess and the observed upper limit of the ra- 
dial velocity impose the upper limit of the mass loss rate 
of several xlO -9 Af®yr _1 (Waters, Marlborough 1994). 

At present, the viscous decretion disk model proposed 
by Lee et al. (1991) is the only model that naturally yields 
near Keplerian disks around Be stars. In this model, 
matter supplied from the equatorial surface of the star 
gradually drifts outward by the viscous torque and forms 
the disk. Assuming Keplerian rotation and neglecting 
the advective term in the equation of motion, Lee et al. 
(1991) obtained steady and thermally-stable structure of 
viscous decretion disks around Be stars. 

Later, using a 3D Smoothed Particle Hydrodynamics 
code, Kroll and Hanuschik (1997) studied the evolution 
of the gas ejected temporarily from a Be star to model 
the transient disk formation/decay process around /i Cen. 
They found that the gaseous particles gradually expand 



and form a near Keplerian disk in the viscous timescale. 

These works show that the viscous torque is likely 
the agent that makes a disk near Keplerian. However, 
since the specific angular momentum in the near Kep- 
lerian disk increases with radius, Be disks are unlikely 
near Keplerian far from the star. A general belief is that 
the disk which is Keplerian near the star becomes an- 
gular momentum conserving far from the star, although 
no theoretical or observational study has confirmed this 
transitional feature of Be disks. 

The purpose of the present paper is to investigate 
whether Be disks have this transitional feature, based on 
the viscous decretion disk scenario. In section 2, we de- 
scribe the disk model and formulate the eigenvalue prob- 
lem for obtaining the radial structure of viscous decretion 
disks. We present the transonic solutions and discuss the 
topology of the sonic points in section 3. In section 4, 
we draw our conclusions and discuss several issues which 
remain unsolved. 

2. Basic Equations for Viscous Decretion Disks 

We assume that the circumstellar disk of a Be star is 
steady, geometrically thin, and symmetric about the ro- 
tational axis and the equatorial plane. For simplicity, we 
adopt the Shakura-Sunyaev's prescription for the viscous 
stress. Since the deviation from a point mass potential 
due to the rotational distortion of the star is small in 
general, we neglect the quadrupole contribution to the 
gravitational potential. We use a cylindrical coordinate 
system (r, <j>, z) to describe the disk. 

Millar and Marlborough (1998, 1999) computed the 
temperature distribution within lOOi? around the BO star 
7 Cas and the B8-9 star 1 Del by balancing at each po- 
sition the rates of energy gain and energy loss, where R 
is the stellar radius. They found that the disk temper- 
ature is roughly constant in the radial direction within 
100i?. Since solving the hydrodynamical equations to- 
gether with the energy equation is a formidable task, we 
assume for simplicity that the disk is isothermal. We be- 
lieve that this assumption is a reasonable one as a first 
step to model the disk. 

Under these assumptions, the vertically-integrated 
equations describing the mass, momentum, and angu- 
lar momentum conservation in the disk around the star 
of mass M, together with the equation of state and the 
Shakura-Sunyaev viscosity prescription, are 
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where M is the decretion rate, E, W, and t rt f, are the 
vertically-integrated density, pressure, and r-<p compo- 
nent of the viscous stress, respectively, V r and V$ are 
the radial and azimuthal components of the vertically- 
averaged velocity, respectively, F Ta d is the vertically- 
averaged radiative force per unit mass, c s is the isother- 
mal sound speed, and a is the viscosity parameter. The 
last term of the left-hand-side of equation (|^) is the cor- 
rection for the decrease of radial component of the effec- 
tive gravity away from the equator (e.g., Matsumoto et 
al. 1984). 

Since the radial flow in the inner part of viscous de- 
cretion disks is considered to be highly subsonic (Lee 
et al. 1991), the radiative acceleration would arise not 
from the optically-thick strong lines but from optically- 
thin weak lines (and from the optically-thin continuum) . 
Since almost nothing is known about the form of the ra- 
diative force due to an ensemble of optically-thin lines, 
we adopt in this paper the parametric form used by Chen 
and Marlborough (1994): 



GMT GAf(l-r) ,r 
rad ~ r 2 + r 2 V \R 



(6) 



where r\ and e are parameters which characterize the force 
due to the ensemble of optically-thin lines. In the above 
expression, T is the Eddington factor that accounts for 
the reduction in the effective gravity due to electron scat- 
tering. In the rest of this paper, we neglect the Eddington 
factor T, because for typical values for Be stars, T is as 
small as - 0.03 for a BOV star and - 0.003 for a B5V 
star. The radiative force is then written as 
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Eliminating W, E, and t r <f, from equations (d)-(H), we 
obtain the following equations which describe the steady 
viscous flow: 
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I = rVrj, is the specific angular momentum, and r s and £ s 
are the radius of the sonic point and the specific angular 
momentum at the sonic point, respectively. It is impor- 
tant to note that the angular momentum distribution of 
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the flow is not given a priori but obtained, with the ra- 
dial velocity distribution, as the solution of equations (j^) 
and (^|). In these equations, r s and l s are the eigenval- 
ues. Note also that, from equation (||), which indicates 
£ + ac 2 r /V r is constant, V r has to approximately linearly 
increase with r in order for the disk to be near Keplerian. 

We solve equations (0) and (0) by the shooting method 
with the starting point at r = R. Each trial integration 
needs two starting values, one of which is freely specifi- 
able. We specify V r as a freely specifiable starting value 
and choose £ so as to be consistent with these equations 
in the very vicinity of r = R. 

3. Viscous Transonic Solutions 

Wind equations (||) and (||) include five independent 
parameters which characterize the solution. They are the 
viscosity parameter a, the isothermal sound speed c s , the 
star's break-up velocity (GM / R) 1 / 2 , and the radiative 
parameters rj and e. Since, in our model, the disk is 
assumed to be isothermal, the structure of the viscous 
transonic outflow is independent of the mass decretion 
rate M [see section 3 of Abramowicz and Kato (1989)]. 
Solving equations (||) and (|J) for a wide range of these 
parameters, we found that a transonic solution exists for 
any combination of parameter values. 

In the first part of this section, we discuss the struc- 
ture of viscous transonic decretion disks. In the second 
part, we discuss the topology of the sonic point. Since 
we found that the decretion disk structure is not sen- 
sitive to the spectral type of the central star, we will 
present only examples in which the central star is a BO 
main-sequence star with M = 17.8 M®, R = 7.41 and 
T cff = 2.80- 10 4 K (Allen 1973), where T eS is the effective 
temperature. 

3.1. Structure of the Viscous Transonic Outflow 

Figure 1 shows the structure of the viscous transonic 
decretion disk for = ^T c g and (77, e) = (0,0), where 
Td is the disk temperature. In figure 1, solid, dashed, 
and dash-dotted lines denote V r /c s , ^/(GM/JJ) 1 / 2 , and 
£/£(i?), respectively. For comparison purpose, transonic 
solutions for a = 1 (thick lines), a = 0.1 (lines with 
intermediate thickness), and a = 0.01 (thin lines) are 
shown. 

From figure 1, we find that the sonic point is located 
far from the star and the outflow is highly subsonic for 
r <C 10 2 R. We also find that physical quantities vary 
approximately as V r oc r, £ oc r~ 2 (i.e., the local den- 
sity oc r~ 7 / 2 in the isothermal disk), and V<f, oc r -1 / 2 
in the inner subsonic part. For smaller a, the slopes of 
V r and E become steeper in the outer subsonic part. In 
the outer subsonic part and in the supersonic part, 
approximately decreases as r -1 . It is important to note 
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Fig. 1.. Structure of the viscous transonic decretion disk for 
T d = \T cS and (rj, e) = (0,0). The central star is a B0 
main-sequence star. Solid, dashed, and dash-dotted lines de- 
note V r /c a , ^/(GM/fl) 1 / 2 , and S/S(K), respectively. Thick 
lines, lines with intermediate thickness, and thin lines are for 
o=l, 0.1, and 0.01, respectively. 



that the transonic decretion disk has the transitional fea- 
ture mentioned in section 1: It is near Keplerian in the 
inner part, while it is angular momentum conserving in 
the outer part. As shown below, these features of viscous 
transonic decretion disks hold irrespective of parameter 
values. 

As mentioned above, there are five parameters that 
characterize the current problem. In figure 1 we have 
studied how the viscosity parameter affects the transonic 
solution with the other four parameters being fixed. In 
figure 2, we show how transonic solutions depend on (a) 
the disk temperature Td/T e g and (b) the radiative pa- 
rameters rj and e. From figure 2, we note that the char- 
acteristics of the transonic solutions are always the same 
as those shown in figure 1. We also note that the sonic ra- 
dius is smaller for higher disk temperature and/or larger 
radiative force. 

It is important to note that it is basically the pres- 
sure force that accelerates the flow up to a supersonic 
speed. The acceleration by the pressure force does not 
work effectively in the region where it is much weaker 
than the effective gravity. This is why the sonic point is 
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Fig. 2.. Parameter dependence of the viscous transonic solutions. The value of a is fixed to be 0.1. The central star is a BO main-sequence 
star, (a) The temperature dependence: Thick lines, lines with intermediate thickness, and thin lines are for Td/T c ff = 1, 1/2, and 
1/10, respectively. No radiative force is included, (b) The dependence on the radiative parameters: Thick lines and lines with 
intermediate thickness are for e) = (0.5,0.1) and (r), e) = (0.1,0.1), respectively. Thin lines are for no radiative force. The disk 
temperature is fixed to be T^/T^g = 1/2. Other format of the figure is the same as that of figure 1. 



located far from the star. From equation (0), the sonic 
radius r s is roughly given by r s ~ iGM/q 2 when the 
radiative force is much weaker than gravity. That the 
pressure force plays an important role for the accelera- 
tion suggests that the flow structure can be sensitive to 
the radial distribution of the disk temperature. We have 
assumed, for simplicity, the disk to be isothermal, but 
the disk temperature should be originally determined by 
solving the energy equation. If the disk is not isothermal, 
it will affect the flow structure, in particular in the region 
of r > 100i?, where the sonic point in the isothermal disk 
lies. 

We would like to make a brief comment here on 
the radial velocity of viscous decretion disks. One 
may estimate from equation (^) that the radial veloc- 
ity V r is of the order of a(H /r) 2 ^, which is reduced 
to 0.6a(T d /T , eff )(r/i?) 1 / 2 kms~ 1 for an isothermal, near 
Keplerian disk around a BO-type main-sequence star. 
As shown in figures 1 and 2, however, this estimate is 
wrong in both absolute magnitude and radial depen- 
dence. Equation ([)]) indicates that V r roughly linearly 
increase with r, as discussed at the end of section 2, and 
has to be much smaller than a(H/r) 2 V,p near the star. 



3. 2. Topology of the Sonic Point 

It is well known that the stability of the transonic ac- 
cretion with Shakura-Sunyaev viscosity prescriptions is 
related to the topology of the sonic point (e.g., Abramow- 
icz, Kato 1989, and references therein). For a smaller 
than a critical value, the sonic point is of saddle type and 
the transonic accretion is stable. On the other hand, for 
a greater than the critical value, the sonic point is nodal 
and the transonic accretion becomes unstable. The in- 
stability arises as a result of work done by viscous force 
[see Kato et al. (1988) for detailed analysis]. 

In the current model, both left-hand and right-hand 
sides of equation (0) must vanish at the sonic point, i.e., 



V r (r s ) = c s 



tl = GMi 



1 — r] 



' 2 2 

-czrt- 



(11) 



(12) 



To examine the topology of the sonic point, we evaluate 
the derivative dVv/dr at r = r s . After some manipula- 
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Fig. 3.. Diagram showing the location of the sonic points in the (r s , A s )-plane for (a) a = 1 and (b) a = 0.1. No radiative force is 
included and A s = i s + ac s r s is normalized by (GMR) 1/I2 . The central star is a BO-type main-sequence star. Each contour denotes the 
relation between r B and A s for a constant sound speed. From the bottommost curve upwards, 7d/T cf j = 1, 0.9, 0.8, . . ., 0.3, 0.2, and 
0.1. The filled circle on each contour indicates the position of the sonic point. In each panel, the dashed line separates the spiral-type 
sonic point region (upper side of the line) from the nodal- type sonic point region (lower side of the line), and the dash-dotted line 
separates the nodal-type region (upper side) from the saddle- type region (lower side). In figure 3b the nodal- type region is so narrow 
that it is indistinguishable on the figure. The thin solid line in each panel indicates a boundary in the (r s , A s )-plane, below which 
the circular orbit at the sonic radius is unstable. 



tions, we have 
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dV r 
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Note that equations ( [l3| ) and (14) are quite similar to 
equations (2.15) and (2.16) of Abramowicz and Kato 
(1989) for the transonic accretion. According to the stan- 
dard classification of critical points (e.g., Ferrari et al. 
1985), the sonic point is of saddle type if X < 0, nodal if 
< X < a 2 £ 2 /4r s 4 , and spiral if X > a 2 4 2 /4r s 4 . 

In figure 3, as typical examples, we present the location 
of the sonic points for (a) a = 1 and (b) a — 0.1. In the 
figure, A s = £ s + ac s r s is normalized by (GMR) 1 / 2 . Solid 
contours denote the relation between r s and A s , given 
by equation @, for T d /T cS = 1, 0.9, 0.8, . . ., 0.3, 0.2, 



and 0.1 from the bottom to the top. The filled circle on 
each contour indicates the position of the sonic point . In 
each panel, the dashed line separates the spiral-type sonic 
point region (upper side of the line) from the nodal-type 
sonic point region (lower side of the line), and the dash- 
dotted line separates the nodal-type region (upper side) 
from the saddle- type region (lower side). In figure 3b the 
nodal-type region is so narrow that it is indistinguishable 
on the figure. The thin solid line in each panel indicates 
a boundary in the (r s , A s )-plane, below which the circular 
orbit at the sonic radius is unstable. 

Figure 3 shows that the sonic point for a = 1 is lo- 
cated in the nodal-type region, whereas that for a = 0.1 
is in the saddle-type region. Examining the transonic so- 
lutions for a wide range of parameters, we found that the 
topology of the sonic point depends only on a and there 
is a critical value, a*, such that the sonic point is nodal 
for a > a», while it is of saddle type for a < a*. Note 
the similarity between the transonic decretion and accre- 
tion. In our model, a» lies between 0.9 and 0.95. This 
means that the steady transonic decretion in disks of Be 
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stars exists for a<0.9, whereas the transonic decretion 
cannot be steady for a > 0.95. 

4. Summary and Discussion 

We have examined the characteristics of the outflow 
in viscous decretion disks around Be stars. For sim- 
plicity, we have assumed the disk to be isothermal and 
adopted the Shakura-Sunyaev a-prescription for the vis- 
cous stress. As the radiative force, we have adopted the 
parametric form used by Chen and Marlborough (1994). 

Solving the resulting wind equations for a wide range 
of parameter values, we have found that a transonic so- 
lution exists for any combination of parameters. When 
the radiative force is much smaller than the gravity, it 
is basically the pressure force that accelerates the flow 
up to a supersonic speed. The sonic radius r s is, then, 
roughly given by r s ~ |GM/c 2 . The sonic radius is 
smaller for higher disk temperature. A larger radiative 
force also gives a smaller sonic radius. The outflow is 
highly subsonic in the inner part of the disk. Approxi- 
mately, the radial velocity increases as r and the surface 
density decreases as r -2 in the inner subsonic part. For 
smaller a, slopes of these quantities become steeper in 
the outer subsonic part. Interestingly, the disk is near 
Keplerian (V^ ~ r -1 / 2 ) in the inner subsonic part, while 
it is angular momentum conserving (y ~ r - v ) in the 
outer subsonic part and in the supersonic part. 

We have also found that the topology of the sonic point 
is nodal for a greater than a critical value, a*, whereas 
it is of saddle type for a < a*. In our model, a* lies 
between 0.9 and 0.95 irrespective of values of other pa- 
rameters. From the similarity between the physics of 
accretion and decretion and the study of the transonic 
accretion disks, we expect that the transonic decretion 
disks with a < a* are stable, whereas those with a > a* 
are unstable. 

We have seen that the viscous transonic decretion disks 
have many characteristics in agreement or consistent with 
the observed features. Thus, the current model is satis- 
factory as a first step to obtain a sophisticated model 
of Be disks. We have to admit, however, that there are 
several issues which remain unsolved by the current ver- 
sion of the model. They are the observed high mass-loss 
rate through the disk or short disk-formation timescale, 
a variety of the observed density distribution, and the 
mass supply mechanism from the star. We would like to 
discuss these issues in the rest of this paper. 

First, we consider the mass loss rate \M\ through the 
disk. In principle, the mass loss rate is determined by 
how much torque is exerted by the star at the inner disk 
boundary. Unfortunately, however, we have no satisfac- 
tory model for the star-disk interaction or observational 
estimate of \M\ in quasi-steady disks (|M| in the disk 
formation stage will be discussed later). Therefore, in 



the below we give an indirect estimate of \M\ the present 
model for the steady, viscous transonic solution allows. 

In isothermal decretion, as in isothermal accretion, M 
is not a parameter which characterizes the problem. It 
does not appear in the wind equations (|§|) and (j^) , so the 
absolute value of the surface density £ is not determined 
from the model. Using the observed base density (i.e., 
the equatorial density at the inner disk radius r = R) to 
normalize S, however, we can estimate M by 



|M| 



2 x 10~ 



Poa(R) V r (R) 



10- n gcm- 3 10- 3 c s 



1/2 



x| 5-) M^yr" 1 . 



(15) 



The base density, poo(R), for most stars analyzed from 
Fe II lines (Hanuschik 1986) and the IR excess (Waters 
et al. 1987; Dougherty et al. 1994) ranges from 10 



-li 



to 



10~ 12 gem -3 . Figure 1 shows that the radial velocity 
at the inner disk radius of the transonic decretion disks, 
V r (R), is roughly of the order of 10 _3 ac s , as far as the 
radiative force is much smaller than the gravity. As a re- 
sult, we obtain \M\ of the order of 10 _11 -10 _12 aM g) yr _1 
for steady, viscous transonic disks. 

It should be noted that this value does not increase 
very much even if we use V r (R) derived from the observed 
upper limit of the disk outflow velocity, i.e., a few kms -1 
for r <J 10i?. Since V r approximately linearly increases 
with r in viscous decretion disks, this upper limit of the 
disk outflow, V r (10R) < a few kms^ 1 ~ 0.1c s , suggests 
that V r (R) < 10 _2 c s and r s >10 2 R. Consequently, using 
the upper limit of the observed radial flow, we have \M\ ~ 
10~ 10 — 10 _11 Af®yr _1 from equation (15). Our results, 
together with the observed range of the base density of 
Be disks, suggest that, in steady viscous decretion, the 
mass loss rate through the disk is at most comparable 
with that in the polar wind region. 

It is interesting to note that another constraint on M is 
obtained by estimating the spin-down timescale of the Be 
star. Porter (1998) derived the constraint on a quantity 
equivalent to M, for which Be stars do not show sig- 
nificant spin-down during their main-sequence lifetimes. 
Porter adopted that each Be star loses its angular mo- 
mentum at the rate of \M\£-k(R), where £k(R) is the Ke- 
plerian angular momentum at r = R (here, for the sake of 
convenience, we have adopted the notation different from 
Porter's). In viscous disks, however, this rate is largely 
an underestimate, because it is likely that the Keplerian 
region extends at least up to lOOi? where the specific an- 
gular momentum is one order of magnitude larger than 
£k(R). Thus, we have to use ~ 10\M\l K (R) instead of 
|M|^k(-R) to evaluate the angular-momentum loss rate 
from the star. Then, the constraint for the star not to 
show significant spin-down during their lifetimes becomes 
more stringent than that derived by Porter (1998). 
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Since it is beyond the scope of the paper to present a 
detailed calculation, we roughly estimate the spin-down 
timescale below. For this purpose, we adopt the same 
stellar model as that adopted by Porter, i.e., a rigidly 
rotating polytrope with the index of 3/2. The stellar 
angular momentum is, then, given by ~ 0.2/M^k(-R), 
where / is the ratio of the stellar rotation velocity to 
the break-up velocity, which ranges from 0.4 for early 
Be stars to 0.8 for late Be stars. Then, the spin- 
down timescale is given by ~ 0.2fM£ K (R)/10\M\e K (R) 
~ 0.02/M/|M| - lO^lMl/lO-^Mgjyr-^yr. Since the 
spin-down timescale has to be much longer than the 
main-sequence lifetime, which is ~ 10 7 yr for early Be 
stars and ~ 10 8 yr for late Be stars, we have the con- 
straint on M for Be stars: \M\ < 10 _8 M®yr _1 for 
early Be stars and \M\ <C 10~ 9 M®yr _1 for late Be 
stars. Since the observed disk structure of Be stars 
is insensitive to the spectral type, the latter constraint 
(|M| < 10" 9 M® y^ 1 ) should be taken. 

Observationally, the rate of equatorial mass loss has 
been derived for several stars only in the disk formation 
stage. Since 1977, fi Cen has exhibited only flickering 
emission caused by episodic mass loss events. Each out- 
burst lasts only for 2-5 d and makes a transient, small 
disk which decays in 20-80 d. Hanuschik et al. (1993) 
estimated the mass loss rate for /z Cen during each out- 
burst to be 4 x 10 _9 Mg,yr _1 . Recently, Telting et al. 
(1998) studied the long-term behavior of the Be disk in 
X Per (4U 0352+30), a Be/X-ray binary system, and de- 
rived the mass loss rate of 4 x 10~ 9 M® yr -1 during a disk 
build-up phase which lasted about ly in 1994-1995. 

These values show that Be stars can eject mass at a 
rate much higher than the upper limit predicted by the 
viscous decretion model for steady disks. This suggests 
that the disk formation is a violent and, consequently, 
non-steady process to which the the steady disk model 
cannot be applied. To study such a process, we need nu- 
merical simulations such as those by Kroll and Hanuschik 
(1997). 

The short timescale of disk formation around these 
stars also suggests that the disk-formation process is es- 
sentially non-steady. For steady, viscous transonic disks 
shown in figure 1, the drift time defined by Jp V~ 1 dr is ~ 
30K.(i?)/10- 3 c s yr at r = 10 R and - 80K.(i?)/10~ 3 c s yr 
at the sonic radius. This means, even if V r (R) is as fast 
as 10~ 2 c s , the formation of an r ~ 10R disk takes about 
3yr and that of a transonic disk takes about 8yr. This 
long drift time suggests that only disks which are persis- 
tent for more than lOyr will have the structure similar 
to those studied in this paper. 

Next, we discuss the density distribution in Be 
disks. In isothermal decretion disks, the radial expo- 
nent of density, n, defined by n = —dlupw/dhir = 
—dln(E/H)/dlnr, is a little bit bigger than 3.5. More 
generally, as pointed out by Porter (1999), 2n + 3m > 7 



must hold in order for the decretion to occur, where m 
is the radial exponent of the disk temperature. Since Be 
disks are roughly isothermal except for the very vicinity 
of the star (Millar, Marlborough 1998, 1999), we would 
expect that the observed density gradient index n is dis- 
tributed around n ~ 3.5 with a small dispersion. 

On the other hand, the observed values of n are dis- 
tributed around n ~ 3 — 3.5, which is in agreement with 
the current model, with a dispersion much wider than 
that predicted. Adopting a simplified disk model, in 
which the disk has a constant opening angle and the den- 
sity depends only on radius, and applying the curve of 
growth method by Waters (1986) to the observed far- 
IR excess for 59 Be stars, Waters et al. (1987) found 
n = 2 — 4. Adopting the same disk model as that of 
Waters et al. (1987), Dougherty et al. (1994) obtained 
n = 2 — 5 from the near-IR excess colors for 144 Be stars. 
Moreover, applying the curve of growth method with a 
more realistic disk model to a Be star, Porter (1999) ob- 
tained a value of n quite similar to that obtained by Wa- 
ters et al. (1987). 

Consequently, the viscous decretion disk model is 
roughly in agreements with the observed density distri- 
bution, but it has to explain why the observed density 
distribution ranges widely. In particular, the model has 
to explain why there are Be stars with n < 3, for which the 
decretion does not seem to occur. Porter (1999) pointed 
out that the viscosity dependent on the radius and/or 
the stellar radiation field may help make a decretion disk 
with a shallow density gradient. More detailed study of 
these possibilities is desired. 

Finally, we come to the most important and long- 
standing issue to be solved, that is, the mechanism which 
enables a star to supply mass into the surrounding disk. 
Be stars rotate significantly slowly compared with the 
break-up velocities. Therefore, in order to supply mass, 
stars have to accelerate matter up to break-up velocities. 
Moreover, since most Be stars keep their disks for long 
periods, either the mass supply must be continuous or 
the interval of mass supply events must be much shorter 
than the viscous timescale. Otherwise, the disk would 
decay before the next mass supply event occurs, like the 
disk of fj, Cen. 

Although there is no widely-accepted mechanism for 
such mass supply, non-radial pulsations of stars may 
work. Osaki (1986) suggested that the dissipation of the 
non-radial pulsations gives matter the necessary angular 
momentum. Recently, Rivinius et al. (1997) found that 
in ,u Cen the time interval of mass supply events is equal 
to a beat period of several strongest non-radial pulsation 
(NRP) modes, indicating that the mass supply from the 
star results from the resonance between strongest NRP 
modes. If this mechanism works in other Be stars, the 
frequencies of the strongest NRP modes in most Be stars 
will be distributed so that the beat period becomes much 
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shorter than the viscous timescale. It is highly desirable 
to study this possibility further. 

The author is grateful for H.F. Henrichs, I. Negueru- 
ela, and J.M. Porter for stimulating discussions. He 
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also thanks the Astronomical Institute 'Anton Pan- 
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